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This paper investigates the statistical properties of the sum, S, of an 
n-dimensional Gaussian random vector, N, plus the sum of M vectors, 
Xi, • • • , Xm , having random amplitudes and independent arbitrary 
orientations in n-dimensional space. We derive expressions for the proba- 
bility density function (p.d.f.) and distribution function (d.f.) of S and of 
its length, |S|, as series expansions involving only the moments of 
|X,|, i = 1, • • • , M. In addition, we find the p.d.f. and d.f. of the pro- 
jection of S onto 1 -dimensional space. Our results are generalizations of 
the n = 2-dimensional problem of finding the statistical properties of a 
sum of constant-amplitude sinusoids having independent uniformly 
distributed phase angles plus Gaussian noise. The latter problem has been 
treated by Rice 1 and Esposito and Wilson, 2 but our results can also deal 
with sinusoids having random amplitudes. When n = 3, our findings 
treat, in the presence of a Gaussian vector, the classical problem of "random 
flights" dating back to Rayleigh. Some calculations for the 2- and 3-di- 
mensional problem are presented, and an application to coherent phase- 
shift-keying communications systems is discussed. 

I. INTRODUCTION 

In a number of problems arising in communications systems, in 
multipath phenomena, and in other areas, the determination of the 
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statistical properties of a sum of sinusoids and Gaussian noise is 
important for evaluating system performance. For this reason there 
has been interest in this problem for a number of years. Rice 1 first 
investigated the statistical properties of the sum of a single constant- 
amplitude sinusoid and Gaussian noise. Later, Esposito and Wilson 2 
considered this same problem but with two constant-amplitude sinu- 
soids having independent uniformly distributed phase angles. More 
recently, Rice 3 studied the properties of a sum of M sinusoids and 
Gaussian noise. In this paper, we look at the natural generalization 
of this problem to n-dimensional space; namely, we determine the 
statistical properties of the sum of an n-dimensional Gaussian random 
vector plus the sum of M vectors having random amplitudes and in- 
dependent arbitrary orientations in n-dimensional space. In the special 
case when n = 2, our results are applicable to the type of problems 
considered by Rice and Esposito and Wilson, but they can also deal 
with any number of sinusoids with random amplitudes. When n = 3, 
our findings treat, in the presence of a Gaussian vector, the classical 
problem of "random flights" dating back to Rayleigh. 

In Section II we give a definition of spherically symmetric random 
n-vectors and state a theorem which characterizes the form of such 
vectors in an n-dimensional spherical coordinate system. We consider 
M independent spherically symmetric vectors, Xi, • • •, Xat, and define 
S = £<iiX,. Using our characterization theorem, we show that the 
even moments, £[|S| 2 *], k = 1, 2, • • •, can be easily expressed in 
terms of only the moments of |X,-|, i = 1, • • • , M. Then with the 
normal vector N ~ 91 (0, <r 2 I) independent of the X.'s, we derive in 
Section III the probability density functions (p.d.f.'s) and distribution 
functions (d.f.'s) of S + N and of |S + N| as series expansions in- 
volving the moments of | S | . In addition, we derive the p.d.f . and d.f . 
of the projection of S -f- N onto 1-dimensional space in terms of a 
similar series expansion. When n = 2 and M = 2, we check that our 
results agree with those of Esposito and Wilson for two constant- 
amplitude sinusoids. 

Last, in Section IV we present some calculations for the 2- and 
3-dimensional problems, and discuss some aspects of the computational 
procedure that we use. Certain of these calculations provide results 
for the probability of error of a binary coherent phase-shift-keying 
communications system operating in the presence of M co-channel 
interferers and Gaussian noise. These results extend previously pub- 
lished computations. 46 Additionally, our findings can be used to find 
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the probability of error of this system operating in the presence of 
more general types of interference. 

II. SPHERICAL SYMMETRY 

The generalization of sinusoids with uniformly distributed phase 
angles are "spherically symmetric" vectors defined in the following 
way (see Refs. 6 and 7) : 

Definition: A random n-vector X = (X h • • •, X n ), n ^ 1, is spheri- 
cally symmetric with matrix o if and only if the co variance matrix of X 
exists,* E(X) = 0, and the joint characteristic function of X is of the 
form : f 

*x(u) = £[V"*'] = MKupu')*] (1) 

for some function h on [0, <*> ) and where p is some n X n (symmetric) 
positive definite matrix .*• Actually, h and p are defined only up to posi- 
tive multiplicative factors. However, in this paper we are only con- 
cerned with spherically symmetric vectors with p = I = identity 
matrix. Then h is uniquely determined and 4> x (u) = h( |u |). We denote 
such a spherically symmetric vector by the notation "X is s.s." 

Note that if Xi and X 2 are two independent s.s. vectors, then clearly 
Xi + X 2 is also s.s. 

The probability density function of an s.s. vector X can be found by 
Bochner's theorem. 8 If /i(|u|) is absolutely integrable, then the p.d.f. 
ofXis: 

px(x) -Ml* I), (2) 

where 

g " {r) " ^ttF" 2 ^^ /o" h M Xn ' 2J l»-*M&r)d\ r > 0, n ^ 1. 

Thus, if X is s.s., its p.d.f. is constant over every n-dimensional sphere 
centered about the origin. This vector is precisely what is meant by a 
"random flight" in a higher dimensional space. 

For our purposes, a more suitable characterization of an s.s. vector 
is given by the following theorem proved in Ref. 9. 



* Expected value will be denoted by #(•)• 

'We denote vectors by boldface characters: u = (ui, ■ ■ -, u„). The character u' 
is the transpose of u. The norm of u is denoted |u | = (uu')*. 

1 For n = 1, a spherically symmetric random variable has an even characteristic 
function, *x(u) = A[p*|u|]. 
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Theorem 1: Suppose X = (X h ■ • -, X n ), n ^ 2, is s.s. Then there exists 
a unique set of random variables R ^ 0, $* £ [0, ir], A; = 1, • • •, n — 2, 
G [0, 2tt] /or w/iz'c/i* 

Xj f = R ( 'it sin # fc ) cos <i> 7 1 ^ j ^ n - 2 
I,-i = R ( "n sin $*) cos (3) 

X» = U ( n sin $ k J sin 5, 

and furthermore (R, <S>i, • • -, $„_2, 5) are independent and have respec- 
tive p.d.f.'s: 

p R (r) = 2tt"' 2 [ r ( 1 ) j' rn ~ l 9n (r) r > 

P * = r ( ^|±± ) ^ [r ( ^ )]"' sin— ^ 

g </>* g 7T (4) 

k = 1, •••,w, - 2 

p»(0) = ^ ^ < 2tt 

for the g n (-) of (2). 

Conversely, if (R, *i, • • • , 3>„-2, 0) are independent and have the 
p.d.f.'s given by (4), and X is defined as in (3), Men X is s.s. 

The utility of this theorem lies in the fact that the random variables 
(R, $i, •••, <£n-2, 6) are independent with specified p.d.f.'s. As an 
immediate corollary, we see from (2), (3), and (4) that: 

Corollary 1: Suppose X = (Xi, ■ ■ ■, X„), n ^ 1, is s.s. Then its p.d.f. 
is given by: 

px(x) = (2^)-»r(n/2) P ( '^^| 1 ) - 
Moreover, forj = 1, • • •, n and for all i, 

r(n/2)#[|x| 2 '] imgCIXjlfj 

T[(n/2) + o r(i + i) 

Using Theorem 1 we can prove: 

Theorem 2: Suppose Xi, • ■ •, Xm are independent s.s. n-vectors, n ^ 1. 
Let Sj = Ei=i X„ j = 1, ■■■, M, and define 



* We define n*-i o* = 1. 
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nP = El\S J n, fc = 0, 1, ••., j = l, ...,M 

if* - JC|Xi|«-l »» =0,1, •••, t«l, ...,Jf. 

Pu* 

a B[(2m + l)/2, (w - l)/2] 
C "' 2m ' B», (n - D/2] 

iw/iere B(-,-) is the beta function . Denote 

V f»> = (Mi 0) ,Mi 2) , •••,MJ 2 "") 
vf* = (ff.fT, •••,^ (2 ' n) ) 

and define D„j to be an (;?? + 1) X (»? + 1) matrix ivhose (k, t)th 
element equals 

BUS) *••££"' *"""■ * «*' 

and is if t < k. 

Then, for j = 2, ■ ■ ■, M, and in = 0, 1, ■ • -, 

M j2m) = £ / 2?W \ C„, 2 ,C n , 2>n . 2 ,- ^ v J2»-20 (g) 

i=0\ « / C„ i2 m 

In matrix form this is 

„(2ra) _ „(2m)r) 

so that 

tftt = v [^D ni2 ---D, lJ _ l D„, j (6) 

for j = 2, -..,JJf. 

Proof: By Theorem 1 we have for each X, a corresponding vector in 
spherical coordinate space : 

X, <-» (Ri, *!.,-, •••,* n _ 2 .,- J 0>). 

Since Sy is a sum of independent s.s. vectors, it is also s.s., so there are 
vectors corresponding to it : 

Si«-> {Pi, Si.,-, ••-, 5- 2 .i,^i). 

Note that MJ 2t) = E[F?~\ and ^ 2m) = »[J^"] and that 

#[cos« *!.,-] = £[cos« gi.i] 

r(n/2) f* ... „_ 2 , 

= r(i)ra(» - D] J. C0Btt8m af/a 

c„.2« if i is even ,_^ 

if i is odd. l 
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Denote the components of S, and X ; as follows : 



so that 

Also, 

if i is odd, and 



Sj = (Sl,j, ■ ■ • , Snj) 

Xy = (Xij, ■ ■ ■, X n ,j), 
Sij = Pj cos $\,j 

X\,j = Rj COS *i,y. 

= EtPflELcoa* kul = (8) 



= wTc,. (9) 

and 

tf[Xg] = #[#f Mcos 2 **!,,-] = »T«..» ( 10 ) 

Noting that Si,,_i is independent of Xij since S/_i depends only on 
Xi, • • • , X/_i which are independent of X j} the following equalities 
follow from (8) to (10) : 

«r°c..ta, = row 



t=0 \ * / 

- Ef^)i[aj-i]«[im 
= £ ( 2 2 "' ) »c*-i]ffcw*] 

« /2m\ M j 



(2m -20 „ . o . 
Cn,2m-2i- 



Hence, 



«r- £(^ V-' e --" /^»> 

,=0 \ 41 / C n ,2m 



(2m -20 



The vector equation 

tf (2m) = v (2_m) I)n . j = 2 , • • • , M (11) 

follows immediately from this expression. Since yP m) = v} 2m) , eq. (11) 
implies that tf* = fTD.,1- ■ Dn,i-xDn,j for j = 2,---,M. For 
n ^ 1, we prove (5) directly. 
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In Reference 9 we proved the result : 

Theorem 3:IfX is s.s. and independent of N ~ 31(0, <r 2 I), then the p.d.f. 
ofX + Nis: 

*>*-«(*) == ^V r ( \ ) f" * P.xi W Mz |](„- 2 )/ 2 

Xexp[-~(t; 8 + |z|*)1j (b - 2 ,, 2 (^\«£ 1, (12) 

where /„(■) is the vth order modified Bessel function of the first kind. 
From this theorem we obtain the following : 

Corollary 2: Suppose X is s.s, independent of N~ 31(0, crl), and 
£«-o[(l/2<r 2 )7iG2S[|X| 2 *] < oc (for example, if |X| is a bounded 
random variable). Let Z be the projection of X + N onto 1-dimensional 
space, i.e., Z is (say) the first component of the vector X + N. Then 
for n ^ 1 the p.d.f.' s ofX + N, of |X -f- N| , and of Z are given, respec- 
tively, by: 



. . r(n/2) / 1 | ,A 

fef.M = (2^7^ 6XP ( " 2^ |Z| V 



v f Lp-^3 ( | z | V2Q ( - l/2g*)'ff[ | X 1 "] 

*« T[(n/2)+t] ' {L6) 



PlmlOO = (2^7-, »"~ l exp(- i t>«) 



v f Lf <-^ (p»/2Q ( - 1/2Q fg | X | M] 

x ,tb r[(«/2) + i] "• U4; 



and 



/ s r(/«/2) / 2 2 \ 



v f Lf-*>frV2g 2 )(-i/2q')'-fl[|x|"] ri « 
x «~ r[( w /2)+i] » (1o) 

where Lf a) (-) are the generalized Laguerre polynomials. In addition, 
for n ^ 1 the "distribution functions" of |X + N| and of Z are given, 
respectively, by: 

P r !|X + N|>a)=^r(|^)-(^)' ,/2 exp(- 2 4 2 ) 

w ^ Lg?(a*/2o*H- l/2^)'g[lX| 2 '] 
*£i~ *T[(n/2) + t] (lb) 

SINUSOIDS AND GAUSSIAN NOISE 563 



and 

x. ,„ , * I a \ oT(n/2) / a 2 \ 

Pr \Z > a) = lerfc(^) ~ 2TWp eXp (~ E 1 ) 

X h tT[(n/2)+»] ' {U) 

where T(-,-) is the incomplete gamma function {Ref. 10, p. 337) and 
erfc (•) is the complementary error function. 

Proof: From Ref. 10, p. 242, we have the generating function 

e Wr U2tx) ~ & r<« + t + i) (18) 

With / = v/Jlo*, x = |z|/V27, and a = (n - 2)/2, we substitute (18) 
into (12) to get: 

Px +N (z) = ^V r ( 1) exp [ - i |z |»] / o M dpp m (f) 

1 ~ ^C(n-2)/2] ( | z | 2/ 2(r a) ( _ 1 / 2g 2) .„2i 

X (2 (r 2)(»-»)/« < tb T[(n/2) + i] 



_ JW2)_ /_J_, |2 \ 
" (2 1 r t r 2 )"' 2eXP V 2** 1 l ) 



*£ r[(n/2)+»] "' Uy; 

assuming that the interchange of integration and expectation is valid. 
The second assertion of the corollary follows from Corollary 1 since 

9-.n/2 

Pix + ni(|z|) = j^/2) |z|- 1 Px + n(z). (20) 

To prove (15) we note that Z = X\ + Ni is a 1-dimensional s.s. 
random vector and we apply eq. (13) (with n = 1) and the second part 
of Corollary 1 to obtain the desired result. 

Next, to show (16) we integrate (14) over the interval (a, «) and 
utilize the relationships : 

/V*-p(-£)*-| <*■>■*■(&£) (2D 
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and for igl, 
/>"-'^(-2v) Lf '"- 2,/2: (^)* 



Equation (21) is given in Ref. 10, p. 337, and (22) is proved in the 
appendix. 

Finally, to obtain (17) we integrate (15) and use eq. (22) with n = 1 
and the definition : 

2 f°° 
erfc (x) = —j / exp ( — t 2 )dt. 

7T* J x 

It remains to justify the various interchanges of integration and 
expectation or summation. For example, to validate the interchange 
in (19) it suffices (Ref. 11, pp. 28-29) to show that 

Since (Ref. 12, p. 207) \L\ a) (y) \ ^ e"' 2 r(a + i + l)/t !r(« + 1), 
the expression in (23) is less than or equal to : 

(l/2<r 2 )< prlvl5l ,„„ / |z|'\ 1 r[(n/2) + i] 



- (l/2«r 2 )'- FrlYl2n /z 2 \l 



T(n/2) 



which is finite by hypothesis. 

The utility of this corollary lies in the fact that we can evaluate the 
various p.d.f.'s and d.f.'s knowing only the moments of | X | and not the 
entire distribution of X. 

III. STATISTICAL PROPERTIES OF THE SUM OF INDEPENDENT 
SPHERICALLY SYMMETRIC VECTORS AND GAUSSIAN NOISE 

For simplicity, we combine the results of Corollary 2 and Theorem 
2 into : 

Theorem 4-' Suppose Xi, • • • , X M are independent s.s. n-vectors, n ^ 1, with 
moments vf m) = #[|X,| 2m ], m = 0, 1, ■ ••, and I = 1, ■■■, M, which 
are also independent o/N ~ 91(0, <r 2 I). Lei S = ££=1 X, and assume that 
]Ci°L (1/2o" 2 )'/jSC|S| 2, Ji ! < oo (for example, if the | X,- 1 's are bounded 
random variables). Let Z be the projection of S + N onto 1-dimensional 
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space. Then the following relations for p.d.f.'s and d.f.'s are valid: 

. . r(n/2) / 1 | l2 \ 

Ps + n(z) = s _ aj - s « q i^-- i | 1 |.j 



X .tb r[(n/2) +0 ' (24) 



/' s - n : ) , . ,-^ji v n ~ l exp ^ - ^i v 2 J 



(2<r 2 )" 



X fa r[(n/2)+i] ' (5) 



, , r(n/2) / z* \ 



- Lf-»>(a'/2a«)(-l/2a«)*Mff ) ^ 

x ,tb r[(n/2) + ' (2b) 



Prf|S + N|>a} =f ^r(|^)-(^) n/2 

Xttp V"2?A5 ir[(»/2) + ' (27) 



and 



t> <v s. i 1 r / a \ aT(n/2) 

X rep f ^ "> t ^WMIZ 1/2^^ (28 ) 
XeXP V V/tt »T[(n/2)+f| W 

77»e moments nffl = i?[|S| 2 *J ore determined by the recurrence relations 
«J» = £ (gn c "-» c "-"-» ^ y J«-"» (29) 

fc=0 \ ^/C / C n ,2i 

for j = 2, ■ ■ -, M with 

_ B[(2m + l)/2, (n - l)/2] 
B[i, (n - l)/2] 

or 6?/ Me matrix equation (6). 

We next look at some special cases. 

A. n = 2-dimensional space 

When n = 2, eqs. (24) through (29) reduce in an obvious manner. 

566 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1974 



The incomplete gamma function in (27) equals (Ref. 10, p. 339) : 

P (l.£)- «p(-£)- (30) 

The Laguerre polynomials L, <-i) and Li* 1 in (26) and (28) can be ex- 
pressed in terms of Hermite polynomials //,(■) using the relations 
(Ref. 10, p. 240) : 

UHx) m <- 1 )'^' ( " ) (31a) 



and 



«-w = ( -=W^- («b) 



For example, eq. (28) can be rewritten as: 
I. r| *>a>=lerfc(^) + Iexp(-£) 



1=1 

We also check that 



v ^ F 2 ,-i(a/V2^)(l/2^)^ ) , Q9 v 



so that 



(2i \ c„, 2 fcC n .2i-2ft _ / i ' \ 
2k ) Cn.ti ~\k) 

-i(S) 

'— ibioG) 2 G) % > 



,(2t)„(2i-2t) 
j/ 2 



k 2 ,{ 2 «^ 2 '- 2t >, 



(2t)(2(-2k) (21-21) 
1 H l> 3 , 

1=0 h=0 \* / \^/ 

and so forth. 

Consider the type of problem investigated by Rice 1,3 and Esposito 
and Wilson 2 in determining the p.d.f.'s of the envelope and instan- 
taneous value of 

u 

2(0 = £ A* cos (w k t + B k ) + n{t), 

where each A k ^ is independent of 0* and 0* is uniformly distributed 
on [0, 27r). Assume that the pairs ( (A k , #*)} are independent of each 
other and of »(•)• Suppose n(t) is the result of the passage of zero- 
mean white stationary Gaussian noise through a bandpass symmetrical 
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filter. Then n(t) can be written as (Ref. 13, pp. 142-148) : 

n(t) = »i(<) cos w t — nt{t) sin wj, 

where ni(t) and w 2 (0 are zero-mean independent stationary low-pass 
Gaussian processes with 

a* = E[n(t)J = El ni (t)J = E[n 2 {t)J. 

Let 6 k (t) «- (wt — w )t + 0* and thus: 

2(0 = Z 4* cos [ua + 0*(O] + »(0 

= Y. A k cos [w £ + 0*(O] + ni(0 cos iv Jl — n 2 (t) sin w t 

A = l 
= £ -4* cos &k{t) + Wi(<) COS «U 

— £ -4* sin Q k (t) + n 2 (0 sin wj, 
= A(t) cos [w t + 7(0]. 
At any time t , let 0* = O k (t ), wi = /ii(< ), and n 2 = n 2 (< P ). Put 
X A = C4*cos0fc, A k smd k ), k = l,---,M, 

and 

N = (m, n 2 ). 

Then Z*^i X* + N is s.s., so by Theorem 1 it has the form (r cos V, 
r sin SP), where r ^ is independent of ^ and SF is uniformly dis- 
tributed on £0, 2x). It follows that 

z(t ) = T cos SF cos Woto — T sin ^ sin wj < o 
= r cos (w t + ^) ; 

that is, r = A(i ) and \F = 7(^0)- Hence, at any time t , A(t ) and 
7 (£„) are independent and y(t ) is uniformly distributed on [0, 2tr). 
Moreover, the p.d.f. of the "envelope" A(t ) is the p.d.f. of 
r ■ I H"=i X* + N| which can be determined from (25). Thus we 
can find the p.d.f. of the envelope of the sum of Gaussian noise plus 
any number of sinusoids with random amplitudes and independent 
uniformly distributed phase angles. The case considered by Esposito 
and Wilson 2 was that of M = 2, A\ = a = constant and A 2 = b 
= constant, in which case 

„p») = .ErjXil 2 "*] = a 2m , 
A 2m) = #[|X 2 | 2 "»] = 6 2 ", 
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and 

— iC) 

The envelope p.d.f. is then: 



r - y i ; ) a ^b^- ik . 



pr(v) =^ 2 exp[-t;V2«r 2 ] 



Li(v 2 /2<r 2 )(- \/2(?y 



x z 



« =o i> I 



L?o(i) w "]' 



which agrees with the result in Esposito and Wilson [Ref. 2, eq. (12)]. 
This expression was obtained earlier by Goldman. 14 

To find the p.d.f. of z(t) at some time instant t , note that : 

z{t ) = A (t ) cos [w t + y(t )~] 

= r cos (w t + &). 

Since * is uniformly distributed on [0, 2ir) and the cosine function 
has period 2ir, the p.d.f. of r cos (w t + ^) is the same as that of 
r cos *. Recall that (r cos *, r sin *) = £f =1 X* + N. Thus, r cos ¥ 
is the first component of the 2-dimensional vector Y,%=i X fc + N and, 
from eq. (26), its p.d.f. is: 

PzU )(zi) = Prco 8 *(zi) 

(- l/2<j*Yntf> 



= (2^p eXP ("2b 2f )lo 



,2 



X «-•'(&)■ (33) 
In Esposito and Wilson's example, this becomes 

(- l/2«*)' 



^ W " (2^ eXp (-27^ ? )fi 



xx * H, (»)[i(i) ,rfVWi ]' 



which agrees with their eq. (29). 

We also check that the d.f. in (27) is the same as that obtained in 
eq. (18) of Ref. 2, when we use the fact that (Ref. 10, p. 241) : 

xLiMx) =ilL?2 1 (x)-Ll»(xn 

Finally, consider a binary coherent phase-shift-keying communica- 
tions system operating in the presence of Gaussian noise and M 
co-channel interferers modeled by a sum of constant amplitude sinu- 
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soids with independent, uniformly distributed phase angles, 6 h • • • , 6m. 
(Details of this model and system may be found in Refs. 4 and 5.) 
The probability of error in such a system is: 4 - 6 

p e = Pr I £ bi cos0« + N x > a\ , 

where iVi ~ 91(0, o- 2 ) and a is the amplitude of the transmitted (de- 
sired) signal. This probability of error is given by the expression in 
(32) and agrees with the result found in Refs. 4 and 5. However, eq. 
(32) can also be used to find the probability of error in this system for 
a more general class of co-channel interference consisting of a sum of 
uniformly phased sinusoids having also independent random ampli- 
tudes. 

B. n — 3-dimensional space 

When n = 3, eqs. (24) through (29) reduce in a straightforward way. 
Equations (24) to (26) and (28) can also be expressed in terms of 
Hermite polynomials by employing eq. (31). The incomplete gamma 
function in (27) can be written in terms of tabulated functions by use 
of the relations (Ref. 10, pp. 339-340) : 

T(c + 1, x) = cT(c, x) + x c e x 
and 

r(i x) = ir*erf (a*), 

where erf(-) is the error function. 

The recurrence relation for the moments becomes, for M = 2, 

and so on for higher values of M. 

IV. SOME COMPUTATIONS 

The form of the expressions in (24) through (28) is quite similar, 
and so the computer programs used for their evaluation were only 
slight modifications of one basic (Fortran IV) program. Different 
values of n could also be treated easily. The basic program required 
computation of a sum of the form : 

f L t <°>(z)(-l/2a*)Vff> 

h r[(n/2)+t] • {6V 

where a; is a variable. 
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2i-2t) 




Fig. 1— Plots of p.d.f. of |S + N | for a 2 = 1, M = 3, n = 2, and different sets of 
vector lengths (bi, bi, b%). 



In one part of our program, the moments nffl were determined from 
eq. (29): 

n} 2i) = t C n (i, fc)MJ 2 -V- 2l) , 3 - 2, • • •, M t (35) 

*=0 



where 



r, /■ ,n A / 2l \c„,2*C„.2i- 



(36) 



Using the definition of c„.2m and properties of the beta function, we can 
show that the coefficients C n (i, ft) are equal to : 



C n (i, ft) 



r(t + l)r(n/2)r[(n/2) + 



r(ft + i)r(i - ft + i)r[(n/2) + ft]r[(n/2) + % - ft] 



(37) 



To efficiently compute these coefficients and to eliminate "overflow" 
problems, we utilized the simple recurrence relation 

n u u\ (t - fe + l)[(n/2) + i - ft] r , . , , x , > , /oon 
C,(t, ft) = ^(n/2) + ft - 1] C " (t ' * 1} ' k ~ h (38) 
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Fig. 2— Plots of p.d.f. of |S + N | for <r J = 1, M = 4, n = 2, and different sets of 
vector lengths (6i, 62, b«, 64). 




Fig. 3— Plots of p.d.f. of |S + N | for a 1 = 1, M = 6, n = 2, and vector lengths 
all equal to 6. 
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Fig. 4— Plots of p.d.f. of |S + N | for * 8 = 1, M = 3, n = 3, and different sets of 
vector lengths (bi, 6j, 6i). 




Fig. 5— Plots of p.d.f. of [S + N | for <r 2 = 1, M = 4, n = 3, and different sets of 
vector lengths (6|, 62, 6j, b t ). 
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Fig. 6— Plots of p.d.f. of |S + N | for <r 2 = 1, M = 6, n = 3, and vector lengths 
all equal to b. 
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Fig. 7— Plots of Pr {Z > a) for n = 2, 10 log 10 (aV££. &i) = 6 dB, 6, 
and for various values of M. 
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Fig. 8— Plots of Pr {Z > a) for n = 2, 10 log, (a 1 / £& 6|) = 8 dB, 6, = b u , 

and For various values of M. 



together with C„(i, 0) = 1. To evaluate nfl? from eq. (35), particular 
sets of moments { vj 2ti \ could be read into the program. However, for 
simplicity we chose spherically symmetric vectors having constant 
lengths 6i, • • •, 6jif- 

The second part of the program was concerned with computation of 

w«(-*)'A(3+0- 

In order to avoid "overflow" difficulties, we actually computed 

tf»WX'/rg + i) with x = - ± ( f 6 ( ) ! . 
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6 8 10 12 14 16 

10LOGi (a 2 /2o-2)INdB 

Fig. 9— Plots of Pr {Z > a) for n = 2, 101og l0 (rf/EJl, &!) = 10 dB, 6, = b u , 

and for various values of M. 

To do this efficiently we used the iterative identity : 

I4*(x)\* _ (2i + a- 1-x) t 

r[(n/2) + i] tT[(n/2) + »] ^" 1 W 

[which follows from the Laguerre polynomial recurrence relation 
(Ref. 10, p. 241)], together with the fact that U a) {x) = 1 and 
L[ a) (x) = a + 1 - x. 

The final part of the program computed the sum : 

y Li " fo) X> a go (40) 

where A#° = m#7(££i &<)" [The factor l/(E£i &0" was built into 
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the computation of eq. (35) in order to find Ai?° J A convergence check 
was provided to end the summation after additional terms did not 
change any significant digits. Although the program was written to 
handle up to 200 terms in the sum, many computations required less 
than 50 terms. As Esposito and Wilson 2 also noted, for certain values 
of x, a 2 , and (&,}, the terms in (40) alternate in sign and have magni- 
tudes of the order 10 16 . For these cases, precision and convergence could 
not be guaranteed. The typical CPU time required to compute eq. (40) 
for 200 values of x was about 10 to 20 seconds in double precision 
arithmetic on the IBM 370/165 system. 

Some representative results of these computations are shown in 
Figs. 1 to 12. Figures 1 to 6 are plots of Pis+ni(«) as a function of v 
for a 2 = 1, for various values of n and M, and for s.s. vectors having 
constant lengths 6i, • • •, 6 at. Curves for n = M = 2 were given in 
Ref. 2. Figures 7 to 12 are plots of Pr [Z > a\ versus the quantity 
lOlogio (a 2 /2<r 2 ) for fixed values of the quantity 10 logw (a*/Y£-i & 2 <) 
and for various values of n and M . In these curves, for simplicity, we 
took 6i = 6 2 = ■ ■ • = bit. As we discussed in the last section, the plots 
in Figs. 7 to 12 represent the probability of error of a binary coherent 
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12 14 16 

10LOG 10 (a2/2o-2)INdB 



20 



Fig. 10— Plots of Pr {Z > a) for n = 3, 10 logio (a 2 /L^i 6?) = 6 dB, b, = • • • = b u , 
and for various values of M . 
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Fie. 11— Plots of Pr {Z > a) forn = 3, 101og 10 (flVE& &') = 8 dB, 6j = 
and for various values of M . 

phase-shift-keying system versus signal-to-noise ratio, 

SNR = 10 logio (a 2 /2<7*) (dB), 

for fixed values of signal-to-interference ratio, 



?M, 



SIR = 10 log™ 



(V!/<) 



(dB). 



These results extend those previously found in Refs. 4 and 5 to larger 
values of SNR and smaller values of SIR. 

V. CONCLUSION 

In this paper we presented expressions for the p.d.f. of a sum of 
spherically symmetric random vectors plus a Gaussian vector in n-di- 
mensional space. We also found expressions for the p.d.f. and d.f. of 
the length of this sum and of the projection of this sum onto 1-dimen- 
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10LOG 10 U2/20-2 ) INdB 

Fie. 12— Plots of Pr {Z>a) for n = 3, 10 logio (oVESi bf) = l0 dB, 6i> 
and tor various values of M. 



>b u , 



sional space. All of these expressions were series expansions involving 
only the moments of the length of the sum of the s.s. vectors. These 
moments could be found from recurrence relations also derived in the 
paper. Some computations of the p.d.f.'s and d.f.'s were presented for 
the 2- and 3-dimensional cases, and an application to a communications 
system was discussed. However, as pointed out earlier in Refs. 2 and 3, 
there are sometimes difficulties in evaluating these p.d.f.'s and d.f.'s 
for certain parameter values, even for the case of s.s. vectors having 
constant lengths. 

APPENDIX 

To prove eq. (22) we use the fact (Ref. 10, p. 241) that, for i £ 1, 

j i le- t t a+i L\'L\ l) {t)'] = ie- l t a U a) {t). 
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Integrating this expression over the interval (a, <x> ) yields 

- e- a a a+l L\l\ l) {a) = if" rHrtj* (t)dt. (41) 

Equation (22) follows from (41) after a simple change of variables. 
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